Virial theorems for trapped cold atoms 
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We present a general virial theorem for quantum particles with arbitrary zero-range or finite-range 
interactions in an arbitrary external potential. We deduce virial theorems for several situations rel- 
evant to trapped cold atoms: zero-range interactions with and without Efimov effect, hard spheres, 
narrow Feshbach resonances, and finite-range interactions. If the scattering length a is varied adia- 
batically in the BEC-BCS crossover, we find that the trapping potential energy as a function of I/a 
has an inflexion point at unitarity. 
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In quantum mechanics, zero-range interactions can 
be expressed as boundary conditions on the many-body 
wavefunction in the limit of vanishing interparticle dis- 
tance These boundary conditions define the do- 
main of the Hamiltonian, i. e. the set of wavefunc- 
tions on which the Hamiltonian is allowed to act. The 
Hamiltonian of a zero-range model differs from the non- 
interacting Hamiltonian only by its domain. In 31?, the 
zero-range model has a long history in nuclear physics 
going back to the work of Wigner, Bethe and Peierls on 
the 2-nucleon problem [§]. 

Zero-range interactions provide an accurate descrip- 
tion of cold atom experiments 0, 0, S|- In particular, 
two-component fermionic atoms in 31? at a broad Fesh- 
bach resonance are well described by zero-range interac- 
tions of scattering length a — 00. This so-called unitary 
limit is completely universal, e. g. the superfluid transi- 
tion temperature is a universal number times the Fermi 
energy flS!- 

A new ingredient in cold atomic systems with respect 
to nuclear physics is the external trapping potential. For 
the unitary Fermi gas in a harmonic trap, the virial the- 
orem 

E = 2E tl . (1) 

was recently shown theoretically and experimentally 0, 
[ToL [TTL EH- Here E is the total energy and E tI is the 
trapping potential energy. 

On the other hand, the traditional virial theorem does 
not concern zero-range interactions, but more usual in- 
teractions described by a potential energy U{f\, . . . ,rjsf), 
where the domain is simply a set of smooth functions. It 
states that the kinetic energy T is one half of the virial: 



(T) 




V ri C/ 



(2) 



for any eigenstate; implying (T) = n/2 (U) if U is a ho- 
mogeneous function of degree n. This theorem is as old as 
many-particle quantum mechanics [HI], and is used e. 
to understand the properties of many-electron atoms [li 



In this paper, we present a general virial theorem for a 
Hamiltonian with an arbitrary domain. In the particular 
case where the domain does not depend on any length 
scale, we recover the virial theorem for the unitary gas 
Eq. |T]) and the traditional virial theorem Eq. ([2]). By 
considering the case of a more general domain, we find 
new virial theorems for several interactions relevant to 
cold atoms: zero-range interactions of arbitrary scatter- 
ing length with or without Efimov effect, hard spheres, 
narrow Feshbach resonances, and finite-range interac- 
tions. Our theorems hold for any trapping potential, in 
any space dimension. They are valid not only for each 
eigenstate, but also at thermal equilibrium provided the 
entropy S is kept constant. For zero-range interactions 
without Efimov effect, the virial theorem implies that for 
any 5, the function _E tr (l/a, S) has an inflexion point at 
the unitary limit l/a = 0. 

General Virial Theorem. Let us consider a quantum 
problem of N particles, with arbitrary statistics and dis- 
persion relations. The position ri of particle i is a vector 
of arbitrary dimension, with continuous or discrete coor- 
dinates. We consider a general Hamiltonian 



H = H' + C/(n 



,r N 



(3) 



where 



H' and its domain depend on p parameters l\, 



1 lp 

which have the dimension of a length, on h, and on some 
arbitrary fixed mass m 

• J7(ri, . . . , r\) is an arbitrary function, which is suffi- 
ciently regular so that the domains of H and H' coincide. 
Then, as shown below: 



E= (u +]^v i -V ri U 



(4) 

q=l H 



for any stationary state of energy E, the partial deriva- 
tives dE/dl q being taken for a fixed function U. 

To derive the above theorem, we use dimensional anal- 
ysis to rewrite U as 

, h 2 X 2 



-/(Ar x , . . . , Ar N ) 



(5) 
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where A has the dimension of the inverse of a length, and 
/ is dimensionless function. The theorem then follows 
from the following two relations: 



dE 



N 



dE 

'ax 



9=1 



(6) 
(7) 



Here the partial derivatives with respect to A are taken 
for a fixed function / and for fixed l\, . . . , l p . 
Eq. © follows from the Hcllmann-Feynman theorem (l5j 
and from Eq. ([5]). The Hcllmann-Feynman theorem holds 
if the derivative d \tp) /OX of the considered eigenstate be- 
longs to the domain of H . We expect this to be true in 
all situations considered in this paper. 
Eq. {7J follows from the fact that, by dimensional analy- 
sis, the energy writes 



E(h,...,l p ,[U}) 



2 \2 



ft* A 



■F(Xh,...,Xl p ,[f]) (8) 



where F is a dimensionless functional. 

The traditional virial theorem Eq. ([2]) is recovered by 
applying the general virial theorem to the case where: 
• The operator H' in Eq. |3]) reduces to the kinetic energy 



Eft » 
- — A r . 
2m, ' 



(9) 



rrii being the mass of particle i; 

• The domain is simply a set a wavefunctions which are 
smooth when particles approach each other. 
Since this domain does not depend on any length scale, 
the second term on the right-hand-side of Eq. ^ van- 
ishes, and the result Eq. ([2]) follows. 
Virial theorems for trapped cold atoms. In what 
follows we restrict to the experimentally relevant case 
where U is a sum of trapping potentials: 

JV 



U(r 1 ,...,r N ) = ^2u i (r i ), (10) 

i=l 

and we rewrite the general virial theorem Eq. (j4]) as: 

1 p dF 

E = 2E tI -±J2 1 ^ ( U ) 

where 

Btrsi^/^rO + iri-VUiCn)). (12) 



2^'dl 



N 



If each Ui is a harmonic trap, then D tr reduces to the 



trapping potential energy: E tt = Y^iLi ( u i( r i)) = E "- 



A. Zero-range interactions. We now assume that 
each pair of particles either interacts via a zero-range 
interaction of scattering length a, or does not interact 
Zero-range interactions are well-known in ID 

2D ui] and 3D @, m, in, s mi m m. 

A.l Universal states. We call universal state a sta- 
tionary state of the zero-range model which depends only 
on the scattering length. All eigenstates are believed to 
be universal in ID and 2D ([if 



171 124| and references 



therein) and in 3D for fermions with two components of 



equal masses 1 31 \ A 

m HE S3, EaM 




32, 



^MEI M M M 0, 0, 

33l . 134 l35l . l36j or unequal 
masses with a mass ratio not too far from one 0, |36| |. 
For 3 bosons in 3D there are both non-universal efimo- 
vian states and universal states [22|, [23| . 

In the Hilbert space generated by universal states, the 
domain of the Hamiltonian depends only on the scatter- 
ing length. Thus Eq. (jlip gives for any universal state: 



E 



1 dE 
2E tr --a—, 

2 oa 



or equivalently 



E = 2D, 



1 dE 
2a 3(1/ a)' 



(13) 



(14) 



This result generalizes the virial theorem Eq. (TT]) to an 
arbitrary scattering length, trapping potential and space 
dimension. Thus it also app lies to quantum gases in low 
dimensions ([37, H, 3^ 4Q| and refs. therein). For the 
case of 2-component fermions in 3 dimensions and power- 
law traps, this result is also contained in two recently 
submitted works: it was found independently by S. Tan 
in J4ll| and rederived using a method similar to ours in 

For a = 00 (which is the unitary limit in 3D and the 
non-interacting limit in ID and 2D), Eq. (TTJJ becomes: 



E = 2E t] 



(15) 



This generalisation of Eq. ([T]) to an arbitrary trap was 
obtained by Y. Castin (unpublished), and is also con- 
tained in the recent independent work of J. Thomas [43j ]. 
Of course it also holds for a = (which is the Tonks- 
Girardeau limit in ID and the non-interacting limit 2D 
and 3D) in accordance with Eq. (1 13[) . 

Taking the second derivative of Eq. (fl"4"|) we obtain: 



a 2 D t , 



d(l/a) 



0. 



(16) 



which means that generically the curve D tr (l/a) has an 
inflexion point exactly at the unitary limit 1/a = 0. 
We can also rewrite Eq. ([M]) in an integral form: 



2 D(a 2 )- ai 2 D( ai ) = -4 



l/a 2 



l/oi 



a 3 E tr {a)d(l/a), (17) 
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which is likely to have a better signal-to-noise ratio than 
Eq. (|14[) when applied to experiments or numerics. 
A. 2 Efimovian states. The boundary condition in the 
limit where two particles approach each other is called 
Bethe-Peierls boundary condition (BPbc). For 3 bosonic 
or distinguishable particles, there exists Efimov bound 
states [19(, and the domain of the zero-range model is 
defined not only by the BPbc in the limit where two 
particles approach each other, but also by an additional 
boundary condition in the limit where all three particles 
approach each other. While the BPbc depends on the 
scattering length a, this additional boundary condition 
depends on a 3-body parameter which we call Rt and 
has the dimensions of a length [23L 13]. The resulting 
2-parameter model is known to be self-adjoint and phys- 
ically meaningful for TV = 3 particles [13, IH 0, El B | • 



The case N > 4 is still controversial [451 ]. 

For this model, the general virial theorem Eq. (fTT|) 
gives: 



E = 2E tI 



1 dE BE 
adlljaj ~ t dR t 



For a = 00 this reduces to 



E = 2E U . 



R t dE 
~2dR t 



(18) 



(19) 



We now apply this to the unitary 3-boson problem in 
an isotropic harmonic trap, which is exactly solvable 22, 
23j ■ The spectrum is E = Eqm +S fku where Eqm is the 
energy of the center-of-mass and £ solves: 



argT 



1 + s-S 



= — |s|lni? t + argr(l + s) mod n, 



(20) 

s ~ i ■ 1.00624 being the only solution s G i ■ (0; +00) of 
the equation: scos(s7r/2) — 8/\/3 sin (sn/6) = 0. This 
allows to calculate d£ /dRt, and Eq. (fl"f)|) then gives [46j |: 



Etr 




l+s-£\ 



(21) 



where ijj is the digamma function. But we can also ex- 
press E tT using the wavefunction, which has a simple ex- 
pression in terms of the Whittaker W function [23j ; the 
result agrees with Eq. (|2"Tj) provided that [46j |: 



dx 



We s 

2 ' 2 



(x) 



£ ■ Ira. ip 



1 



£ 



2 7T 



sinh(|a|7r) |r(i±§=£) 



(22) 



Numerical checks confirm this relation. 

B. Hard sphere interactions. Here the domain is 

defined by the condition that the wavefunction vanishes 



if any interparticle distance is smaller than a. Applying 
the general virial theorem with a single length scale gives: 



E = 2E tI -\a^. 

2 oa 



(23) 



Again, it can be useful to rewrite Eq. (f2"3"|) in an integral 
form: 



E(a) = -J / a'E tr {a')da'. 



(24) 



Within the 3D Gross-Pitaevskii theory, adE/da is the 
interaction energy, so that Eq. (|23[) agrees with the virial 
theorem of [4^ ]. 

C. Finite-range interactions. We now consider mod- 
els with two parameters, the scattering length a and a 
range I. Popular examples are the square- well interaction 
potential [28j , separable potentials [23j , and Hubbard-likc 
lattice models where the lattice spacing / plays the role 
of the interaction range @, [H, [5^ . For such 2-parameter 
models the general virial theorem gives: 



E — 2E tI 
and for a = 00: 



1 dE dE 
ad(l/a) ~ ~dl 



I dE 



E - 2Et ^2dV 
Setting E = lim ; ^ E{1), Eq. ([2H]) implies 

E = SE-AE ti . + O{l 2 ) 1 



(25) 



(26) 



(27) 



which can be used to compute numerically Eq. This 
method is simpler than the usual one, where one com- 
putes E for several values of I and extrapolates linearly 

to i =0 minm. 



D. Effective range model and narrow resonances. 

The effective range model has two parameters, the 
scattering length a and the effective range r e . For 
r e < 0, the model describes a narrow Feshbach reso- 
nance 0, |H, H 0, [Hl[. For r e 0', the model has a 
limit cycle described by the zero-range model of Sec. A. 2, 
with R t = Cr e , where the constant C was obtained nu- 
merically [iH and analytically [H(|. The model is ex- 
pected to be hermitian for a modified scalar product, for 
2 particles [H2] and 3 particles (4^|. Thus the Hellmann- 
Feynman theorem can be used and the general virial the- 
orem holds, implying: 



E = 2E*, 



1 dE 
ad(l/a) 



dE 

' dr e 



(28) 



For r e > 0, the effective range model is well-defined if 
r e is treated perturbatively [44|, and Eq. (|28|) holds, in 
agreement with Eq. |25|) . 

At finite temperature. We will show that the above 
results remain true at finite temperature, provided one 



4 



considers adiabatic transformations. For concreteness 
we restrict to zero-range interactions in the universal 
case. We consider that each eigenstate n has an oc- 
cupation probability p n . We set E = J2 n EnPn and 

Let us first recall the reasoning of Tan [4l|, |53j: for 
a finite system, in the limit where a is varied infinitely 
slowly, the adiabatic theorem implies that the p„'s re- 
main constant, so that 



E 



#)• 



d 



d(l/a) 



22E n p n . 



(29) 



Tan concludes that E and E tI can be replaced by their 
average values E and E tr in the virial theorem Eq. (Tj"4")) . 

Alternatively, let us assume that the p n 's are given by 
the canonical distribution p n cx e~ E ^w/i k BT) ^ wnere the 
temperature T varies with a in such a way that the en- 
tropy S = — fcs J2 n Pn l n -P" remains constant. According 
to the principles of thermodynamics, this assumption is 
a good effective description of adiabatic sweep experi- 
ments where a is changed at a rate much smaller than 
thermalisation rates and much larger than heating and 
evaporation rates [g, |3fJ, 35l 154 , l55j . Under this assump- 
tion Eq. ||5SJ) also holds [Hi]. Thus: 



E = 2 Etr + — 



dE 



2a \d(l/a) 



(30) 



This result is physically consistent with Tan's conclusion. 
Moreover it implies: 



d 2 E 



d{l/a 



" '-=0,5 1=0, 



(31) 



a 2 2 E(a 2 , S)- ai 2 E(a u S) = -4 



l/a 2 



l/oi 



a 3 E tr (a,S) d{l/a). 

(32) 

Experimental considerations. Both E and E tr are 
measurable. Indeed, E tr and the trapping potential en- 
ergy E tT can be deduced from an in-situ image of the den- 
sity profile 32,34, 3^, 57 1 , and the released energy E—E tI 
from a time of flight image 33 , H, 58|, [59[ . By measur- 
ing E and E tI , and using the virial theorem Eq . (|3T))) . 
one could deduce the quantity (dE/d(l/a))s [6l|. This 
quantity is also related to the large-momentum tail of the 
momentum distribution [53[ and to the total energy [H3] . 

Moreover, Eqs. (|30|31|32|) can be directly checked by 
measuring E(a) and E tT (a) in an adiabatic sweep exper- 
iment. 
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